Astronomy Sz Astrophysics manuscript no. depaolis 


February 2, 2008 


(DOI: will be inserted by hand later) 





A super massive black hole binary in 3C66B: future 
observational perspectives 



F. De Paolis, G. Ingrosso, A. A. Nucita 



o 
o 

(N 



> 
oo 

On 

t> 
O 

o 
6 

in 
c3 



Dipartimento di Fisica, Universita di Lecce, and INFN, Sezione di Lecce, Via Arnesano, CP 193, F73100 Lecce, 
Italy 

the date of receipt and acceptance should be inserted later 

Abstract. Supermassive black hole binaries (SMBHBs) may exist in the centers of galaxies and active galactic 
nuclei (AGN) and are expected to be fairly common in the Universe as a consequence of merging processes 
between galaxies. The existence of SMBHBs can be probed by looking for double nuclei in galaxy ce nters or, more 
easily, detecting periodic behavior in the observed radio light curves. In a recent paper, Sudou et al. lSudodll2003h 
announced the first direct observation of an SMBHB. Using VLBI observations they found that the unresolved 
radio core of the radio galaxy 3C66B shows a well defined elliptical motion with a period of 1.05 ± 0.03 yrs, 
implying the presence of a couple of massive black holes in the center of the galaxy. In the present paper we 
study the astrophysical implications of the existence of such an SMBHB in 3C66B. In particular we focus on the 
information that can be obtained from detecting a signal periodicity either in the X-ray and/or 7-ray light curves 
as a consequence of the motion of the black holes. These observations could be used to extract further information 
on the physical parameters of the SMBHB and partially solve the system parameter degeneracy. The detection 
of the gravitational wave spectrum emitted by such system may be used to completely determine the physical 
parameters of the binary system. 
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1. Introduction 

It is now well established that galaxies genera lly contain 
massive blac k hole s in their nuclei (Reesll984l Kormendy 
> ■ k Richstone ll995l Richstone et al.[l998) and since mcrg- 
k> , ing ev ents a re expected to be rather common (see e.g. 
"2 ■ White [l993) supermassive black hole binaries (hereinafter 
. 53 ; SMBHBs) may be fo und frequently in galactic cores 
(Begelman et al. ll980|h 

Up to now, several astrophysical phenomena have been 
attributed to SMBHBs. In fact, their presence influences 
the stars in the galaxy center regions whose mass density 
profile, as obs erved by numerical simulations (Nakano & 
Makino Il99l|h may show a p oc r ' profile typical of 
giant galaxies that have experienced many merging pro- 
cesses in their lifetime. 

Moreover, if the black holes are rotating and have mis- 
aligned spins, plasma beams (aligned with the black hole 
axis) precess with period of 10 3 — 10 4 yr and curved jets 
should be observed. This phenomeno n happ ens in double- 
double radio galaxies fSchoenmakers EoOOl) and in both S 



an d X-sh aped galaxi es (Le ahy , 
a.1. ll9S,4 Wang et al. l200,^ . 



Williams 1984, Parma et 
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Periodic outburst activity, as for example in the quasar 
OJ 2 87 (see e.g. Sillanpaa et al. Il998l Lehto & Valtonen 
1996) and periodic behavior in the radio, optical, X-ray 
and 7-ray light curves of many AGNs are also possibly 
related to the presence of a central massive binary black 
hole creating a jet either aligned along the li ne of sight or 
interacting with an accretion disk (Yu|2002). 

Up to now, the search for X-ray and/or 7-ray vari- 
ability was considered as a method to probe the existence 
of an SMBHB in the center of a galaxy. For example, at 
least three Mkn objects (Mkn 501, Mkn 421 and Mkn 766) 
reveal periodic behavior in the observed signal so that in 
recent studies they have been considered as pos sibly con- 
taini ng SM BHBs (see Rieger & Mann heim koOOL De Paolis 
et al. I2OO2I and De Paolis et al. l2003|h 

However, the best tracer of the existence of SMBHB 
in a galaxy core remains the detection of the Keplerian 
orbital motion of some emission component close to the 
two black holes and, eventually, the observation of eclipse 
episo des as in the case of the AGN OJ 287 (Sadun et al. 
11994 Lehto & Valtonen Il99l Pietila ll998(l . 

Unfortunately, a technical problem is that the typical 
black hole separation (~ 10 16 — 10 cm), for a host galaxy 
located at the distance of ~ 100 Mpc, corresponds to a 
separation of ~ 10 — 100 /Ltarcsec and only recently has 



2 



De Paolis et al.: A SMBH binary in 3C66B: future observational perspectives 



the required accuracy in position measurements become 
available by the phase-referencing very long baseline in- 
terferometry technique at radio frequencies. 

In a recent paper Sudou et al. (|2003|) . using the VLBI 
facility, announced the first detection of a well defined 
Keplcrian orbital motion of the radio emission compo- 
nent in a pair of well separated radio sources, 3C66B and 
3C66A. The natural explanation of this observation is that 
the radio galaxy 3C66B hosts an SMBHB at its center. 

On the as sump tion that the SMBHB orbit is circular, 
Sudou et al. (2003) have estimated, by a best fitting pro- 
cedure, the physical parameters of the system and found 
a mass density value of 2 x 10 15 Mq pc ~ 3 , larger than of 
the massive object in NGC4258 and SgrA*. 

However, radio observations do not make it possible 
to determine all the binary system parameters since the 
results depend on the mass ratio q between the two com- 
ponents, on the eccentricity of the system, and on the 
semi-major axis. 

The aim of this paper is to present a complete analysis 
of the SMBHB in question by relaxing the assumption of 
a circular orbit 1 and indicating some direction for further 
observations which may solve the parameter degeneracy 
for the system in 3C66B. 

In fact, a galaxy nucleus hosting an SMBHB of which 
one of the components is emitting a jet with Lorentz 
factor 7l towards the observer may show periodic light 
curves both at X-r ay and 7-ray wavel engths (see Rieger 
& Mannheim I2OO0I and De Paolis et al. l2002|) . 

Thus, assuming that the radio jet in 3C66B is emitted 
by the less massive black hole (with mass to) which is 
moving in an elliptical orbit (with semi-major axis a given 
by the VLBI radio observations, and generic eccentricity 
e) around its companion (having mass M), we study both 
the expected AT-ray and 7-ray signal periodicity P bs and 
the flux ratio / between maximum and minimum signal 
as a function of the jet Lorentz factor 7l and of the power 
law spectral index a of the photon flux. 

We show how the measurements of / and P bs may 
be used to extract further information about the physical 
parameters of the SMBHB system hosted by the radio 
galaxy 3C66B. 

We also consider the gravitational wave spectrum and 
amplitude emitted by the SMBHB as a function of the or- 
bital parameters of the system and compare the expected 
signal with the capabilities of the next generation space- 
based interferometers LISA and ASTROD. 

The paper is structured as follows: in Sect. 2 we re- 
view the main properties of the radio galaxy 3C66B and 
in particular the observed parameters of the SMBHB at 
its center. In Sect. 3 we discuss the expected periodicity of 
the light curves in the X-ray and/or 7-ray band and show 
how the flux ratio / depends on the system parameters. 
We show how one can use the expected periodicity P bs 



and / to constrain the physical parameters of the binary 
system. We have performed this analysis by assuming that 
the jet inclination angle to the line of sight is i ~ I/7L as 
expected for typical blazars. In addition, we have consid- 
ered also the more realistic case in which i may be esti- 
mated by radio observations. However, in both cases the 
SMBHB parameters may be completely determined only 
if the spectrum of the gravitational waves emitted by the 
system is detected. Finally, our conclusions are given in 
Sect. 4. 

2. An SMBHB in the radio galaxy 3C66B 

Using the very-long-baseline-interferometer, Sudou et al. 
(2003) recently announced the first direct evidence of the 
existence of an SMBHB in a distant galaxy. 

They looked for the Keplerian orbital motion of a ra- 
dio emission component in a pair of well separated (~ 6 
arcmin) radio sources, 3C66B and 3C66A; 3C66B i s a ra- 
dio galaxy at redshift z = 0.0215 (Matthews et al. 1964) 
where as 3C 66A is a more distant BL Lac object (z = 0.44, 
Miller ll97Sl) which was used as the stationary position ref- 
erence to 3C66B. In fact, observ ing both objects at 2.3 
and 8.4 GHz (Sudou et al. l2003j) the radio core position 
of 3C66B at each epoch was measured with respect to 
3C66A. This revealed a time variation of the 3C66B core 
position which is well fitted by an elliptical motion with 
an averaged period esti mated to be 1.05 ± 0.03 yrs (see 
Table [2 and Sudou et al. l2003l for more details). 

Since the radio core is located at the root of the jet 
where the optical depth for synch rotron self-absorption 
is of the order of unity fLobanov Il998h . it is expected 
that the 2.3 GHz core is located at a greater distance 
from the central engine than the 8.4 GHz core. Thus, the 
observation that the 2.3 GHz semi-major axis (243 ± 30 
yuarcsecs) is larger than the 8.4 GHz one (45 ±4 /^arcsecs) 
is a proof that the jet in 3C66B als o suff ers a precessional 
motion. According to Sudou et al. 2003, the most plausi- 
ble explanation for this result is that while the 2.3 GHz 
core motion maps the jet precession, the elliptic motion 
observed in the 8.4 GHz band well describes the orbital 
motion of an SMBHB in 3C66B with a Keplerian period 
P K = 1.05 ±0.03 yrs. 

Since, as stated above, observations at 8.4 GHz are 
believed to reflect the radio core motion of the SMBHB 
more closely than those at 2.3 GHz, assuming that the 
jet is emanating from the le ss massive black hole (of mass 
to < M), Sudou et al.[2003 estimated the semi-major axis 
a of the fitted orbit to be the upper limit of the orbital 
radius of the SMBHB, i.e. 



5.4 x 10 16 (1 +q) cm, 



(1) 



where q = m/M is the mass ratio between the two binary 
components. With this assumption, and from the defini- 
tion of the Keplerian period 



1 It is expected in fact that if SMBHBs come from merging 
events between galaxies they could be in ell iptical orbits with 
eccentricities up to 0.8 - 0.9 (Titchett 119871) . 



Pk = 27T, 



G(M + to) ' 



(2) 
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Table 1. Parameters of the fitted elliptical motion of the 
radio core in 3C 66B monitored in the two radio band s 2.3 
and 8.4 GHz. Data are taken from Sudou et al. 120031 



Radio Bands 


2.3 GHZ 


8.4 GHZ 


Semi-major axis (/^arcsec) 
Axial ratio 
Period (yrs) 


243 ± 30 
0.31± 0.17 
1.10 ± 0.06 


45 ± 4 
0.24 ± 0.14 
1.02 ± 0.04 



the following relation holds 
M~4.6x 10 10 (l + g) 2 Mq 



(3) 



which makes it possible to estimate the masses of the two 
components once the mass ratio q is given. In Fig. ^ the 
two black hole masses m and M of the black holes are 
shown as a function of the mass ratio q. 



io 12 



10'° 



Fig. 1. The masses of the two black holes m (dashed line) 
and M (solid line) are given as a function of the mass ra- 
tio q. Red lines correspond to the orbital maximum semi- 
major axis a ~ 5.4 x 10 16 (l + g) cm at 8.4 GHz while black 
lines correspond to a ~ 2.9 x 10 (1 + q) cm at 2.3 GHz. 



As is clear from Table ^ the estimated orbit semi- 
major axis a is different at different wavelengths. In fact, 
it is a max ~ 5.4 x 10 16 (1 + q) cm at 8.4 GHz and increases 
by a factor 77 ~ 5.4 at 2.3 GHz becoming a max — 2.9 x 
10 17 (1 + <?) cm. Thus, the mass M in the previous equation 
is increased by a factor rf (as is also clearly stated in 
Sudou et al. l2003f) if the semi-major axis is taken to be 
given by the observations at 2.3 GHz. 

Here we would like to emphasize that observations at 
higher frequencies should allow to map more accurately 
the core motion in 3C66B giving a better estimate of the 
orbit semi-major axis which is expected to be smaller than 
a max ~ 5.4 x 10 16 (1 + q) cm. 



Sudou et al. ( 2003) assumed that the binary system is 
characterized by a circular orbit (i.e. e = 0). However, if 
binary black holes originated from galactic mergers, they 
could be in eccentric orbits, and eccentricity value s up 
to 0.8-0.9 are not necessarily extreme 

(Fitchett|l2s3)- Of 

course, due to gravitational wave emission, orbits tend 
to circularize but this happens within a time-scale of the 
same order of magni tude as the merging time-scale (Peters 
Il964l Fitchett ll987h . Therefore, if an SMBHB is found at 
the center of a galaxy, it may happen that the compo- 
nents are still in eccentric orbits. Consequently, the axial 
ratio quoted in Table ^ gives an upper limit to the true 
orbital eccentricity e since the binary orbit is seen pro- 
jected on the plane of the sky. Thus, the possible SMBHB 
in 3C66B may be characterized by an orbit with with an 
eccentricity up to e max ~ 0.96. Relaxing the assumption 
of circular orbits will produce major changes in the esti- 
mation of both the gravitational wave amplitude and the 
coalescing time-scale of the SMBHB possibly hosted by 
3C66B. 



3. Signal periodicity in the X-ray and 7-ray light 
curves and the determination of the 
maximum-to-minimum flux ratio 

Periodic outburst activities and periodic behavior in the 
radio, optical, X-ray and 7-ray light curves of many AGNs 
are believed to be related to the pr esence of an SMBHB 
in the center of the host galaxy (Yu 120021). 

On this basis, Rieger & Mannheim j2000) (but see 
also De Paolis et al. 2002 for a generalization to orbits 
with generic eccentricity values) have recently proposed 
a method to determine the orbital parameters of the bi- 
nary system from the observed quantities, i.e. the signal 
periodicity fobs, the flux ratio / between maximum and 
minimum signal and the power law spectral index a of the 
photon flux. 

According to these models, the periodicity in the flar- 
ing state is assumed to be the consequence of the or- 
bital motion of a relativistic jet in the binary black hole. 
Therefore, the observed signal periodicity has a geometri- 
cal origin due to Doppler-shifted modulation. 

As a consequence of the binary orbital revolution 
around the center of mass, the expected signal period 
Pobs is rel ated t o the Keplerian period Pk by (Rieger & 
Mannheim boOOft 



Pobs = (1 



1 - 



(4) 



where, according to th e stan dard model for the jet bulk 
motion (see e.g. Spada fT999|h the jet outflow velocity v z 
is given by 



o(l ~ lht) 1/2 



(5) 



and i is the inclination angle between the jet axis and 
the line of sight. In the case o f blaz ars, it is usually as- 
sumed that i ~ CRobson Il996t Urry and Padovani 
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1995). However, in the case of radio galaxies this assump- 
tion may not always hold and i (or at least an upper limit) 
can be estimated by direct observations. In the following 
we develop the necessary formalism in general (without 
assumptions for the inclination angle) and then we will 
apply it to the radio galaxy 3C66B by assuming either 
i ~ 1 /tl or the available radio observa tions w hich restrict 
the parameter range (Hardcastle et al. 1996). 

For a typical blazar in which i ~ 1/tl, we show in Fig. 
13 the dependence of the expected light curve periodicity 
on 7 L . 

Thus, for the 3C66B galaxy it is expected that in the 
X-ray and/or 7-ray band the signal periodicity P bs cor- 
responding to Pk — 1-05 yrs is in the range ~ 10 — 150 
days as a function of 7l in the case of i ~ 1/7l- 3C66B 




2 4 6 8 10 



7l 

Fig. 2. The expected light curve periodicity as a function 
of the jet Lorentz factor 7l is shown assuming the relation 
i ~ 1/tl- 



is a large source (several arcminutes in the direction of 
largest angular size) and is characterized by a jet and a 
counterjet with a high degree of asymmetry in the surface 
brightness up to distances of a few kpc from the radio 
nucleus, becoming more symmetric on largest scales. 

The flux asymmetry may be explained by consider- 
ing relativistic beaming, i.e. the whole jet material moves 
at the same bulk velocity (3c and the jet and counterjet 
make the same angle i to the line of sight. In this case, 
the emission is Doppler-beamed towards or away from the 
observer and the ratio R between the jet a nd co unterjet 
radio emission is given by (Hardcastle et al.^996) 

where a r is the spectral index in the radio band. 
Measuring the integrated flux ratios of the jet and coun- 
terjet at different angular distances, it is possible to find 



(3cosi as a function of distance from the nucleus. This 
procedure constrains (3cosi to be in the range 0.25 — 0.6. 
Moreover, values of /3cosi around 0.6 constrain the angle 
to the line of sight of the emerging jet to be i < 53°. For 
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Fig. 3. The expected light curve periodicity as a function 
of 7l is shown for different jet viewing angles i. 



our purposes, leaving the jet inclination angle i as a free 
parameter for the bulk motion velocity (3 one has the ob- 
vious condition (3 m i n < (3 < /3 m ax, with f3 m i n = 0.25/ cosi 
and (3 max =0.6/ cosi. Consequently, the Lorentz factor is 
in the range (1 - P 2 min )' - 5 < 7 L < (1 - ALJ" ' 5 - 

In this way we can evaluate the expected X-ray and/or 
7-ray light curve periodicity as a function of 7l for differ- 
ent jet inclination angles i by using Eq. (@J. In Fig. 13 we 
plot P obs for the angles i = 20°, 30°, 45° and 50°. Note 
that the observed periodicity is always less than one year. 

Let us further assume that the jet is emitted by the 
less massive black hole (of mass m) and that, as usual, 
the jet axis is parallel to the orbital angular momentum 
vector. 2 Moreover, we assume that the non-thermal X- 
ray and/or 7-ray radiation propagates outwards from the 
core along the jet with Lorentz factor 7l. In this way, the 
observed flux modulation due to Dopp ler boosting can be 
written as (Rieger & Mannheim 2000) 

S(u) = 6 3+a S'(v) , (7) 

2 We would like to mention that a generalization of our treat- 
ment to the more general case with arbitrary direction of the 
jet axis with respect to the orbital angular momentum vector 
requires the introduction of a set of quantities (the position 
angles of the jet in the sky plane, the jet precession opening 
half-angle, the inclination angle of the orbital plane, etc.) fol- 
lowing, for example, a simi lar approach to that proposed by 
Abraham & Carrara However, this approach requires 

the knowledge of a number of parameters that are only weakly 
constrained by the observations. A comprehensive treatment 
of this case will be presented elsewhere. 
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where a is the source spectral index 3 and the Doppler 
factor is given by 



^1-( V 2 +Vls (0)2)/ C 2 

1 — (v z cos i + vis (9) sin i)/c 



(8) 



Here, vi s is the component of the less massive black hole 
velocity along the line of sight and 9 the polar angle defin- 
ing the position of m in its orbit. 

As discussed in De Paolis et al. l |2002l) . the maxi- 
mum and minimum values of the velocity vi s are attained 
for 8 — 7r/2 and 9 = 3n/2, respectively corresponding, 
through Eq. ||SJ, to the two possible values S max and S m i n 
of the Doppler factor. With the assumption that the peri- 
odicity in the observed signal is due to the orbital motion 
of the binary black hole, one obtains from Eq. (0) the 
condition 



f 



1/(3+") 



(9) 



where / is the observed maximum to minimum flux ratio. 
Once the value of vi s is known the pr evious equations give 
(for details see De Paolis et al. [2002) 



M 



P. 



1/3 



obs 



(rn + M) 2 / 3 [2tt(1 + z)G]V3 



yl/3+c 



(10) 



1 



Jl/3+a _|_ I 



1 cos i 

c 



2/3 



(1 



=2U/2 



VLBI observations may be used to constrain the 
masses of the two black holes in 3C66B so that the previ- 
ous equation can be solved with respect to / giving 



/ 



with 



k = 



1 



3+a 



(11) 



M 



[2tt(1 + z)G} 1 / 3 sini 



(to + M) 2 /3 



P 



1/3 



obs 



(12) 



(l-^cosi) 



-2/3 



1-— cosi) (1 - e ) 



-1/2 



The percentage flux variation given by the previ- 
ous equation depends on the orbital parameters of the 
SMBHB systems, the jet viewing angle i and the jet 
Lorentz factor 7l. As discussed above, for a typical blazar 
the latter quantities depend on each other (i.e. i ~ 1/7l) 
while for radio galaxies like 3C66B this is not always true. 

In Fig. 0| we show the expected percentage flux vari- 
ation between the maximum and minimum signal as a 
function of the orbital eccentricity, for different values of 
both the Lorentz factor 7l and the power law spectral 
index a. Here we have considered the galaxy as a typical 
blazar, thus setting i ~ 1/tl in ah the previous equations. 



In the following, typical values of the power law index a 
are a ssume d to be in the r ange 1 -2 (for details see Guainazzi, 
et al. ll999l and Boiler et al. l2000l) . 



Since we expect that future radio observations at fre- 
quencies higher than 8.4 GHz will improve the mapping of 
the SMBHB orbit in 3C66B, we consider, as an example, 
the orbit semi-major axis to be a ~ 5.4 x 10 16 (1 + q)jr\ 
cm, being rj ~ 4. 

Note that in this case, using the orbit semi-major axis 
obtained by the 8.4 GHz radio map (i.e. a max ~ 5.4 x 
10 16 (1 + q) cm), Eq. will be satisfied only for small 
Lorentz factor values 7l < 2. In this case, the expected flux 
ratio / is high enough to make 1 — 1// approach unity for 
any value of the eccentricity e (and therefore always imply 
/>100). 



0.7 r 



0.5 : 

0.0 




Fig. 4. The expected percentage flux variation of the X- 
ray and/or 7-ray light curves as a function of the or- 
bit eccentricity. The curves have been obtained assum- 
ing, for the black hole maximum semi-major axis, a ~ 
5.4 x 10 16 (1 + q)/ri cm with r\ ~ 4. Black, blue, red and 
green lines correspond to Lorentz factors 7l of 2, 4, 6 and 
10, respectively. Solid, dotted and dashed lines correspond 
to power law indices a — 1, 1.5 and 2, respectively. Here 
we have also assumed i ~ 1/7l- 



Let us now consider again the most realistic scenario 
according to which the jet inclination angle is i < 53° 

(Hardcastle et al. [l996). In this case, for each possible 
value of the inclination angle i there exists a different set 
of allowed values of the Lorentz factor 7l (see discussion 
above). Hence, using Eq. (|llfl we can evaluate the ex- 
pected flux percentage variation for different jet inclina- 
tion angles i. This is shown in Figs. [S] a) and b) where we 
plot the expected percentage flux for two selected values 
of the jet inclination angle and 7l (i = 50° and 7l = 2 in 
Fig.Et, i = 30° and = 1.3 in Fig.Et). 

As one can see, X-ray and/or 7-ray observations to- 
wards 3C66B may make it possible to detect a possible 
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Fig. 5. The expected percentage flux variation of the X-ray and/or 7-ray light curves as a function of the orbit 
eccentricity. The curves were obtained assuming, for the black hole semi-major axis, a ~ 5.4 x 10 16 (1 + q)/2 cm (black 
lines) and a ~ 5.4 x 10 16 (1 + g)/4 cm (red lines). Solid and dotted lines correspond to power law index a — 1 and 2, 
respectively. Here we are considering a more realistic geometry for the 3C66B galaxy jets. 



periodic behavior in the light curves which can be di- 
rectly used to fix the Lorentz factor 7l (see Fig. EJ. In 
addition, the detection of the signal flux ratio / will al- 
low to determine the orbit eccentricity e of the SMBHB 
possibly hosted by the radio galaxy 3C66B at least if the 
percentage flux variation is not close to unity, and also 
give information about the jet inclination angle i. 

Note also that radio observations of the nuclear region 
of 3C66B have shown a strong flux variability (up to 35 
per cent at high frequencies) on a time scale of months 
(Hardcastle et al. ll996|) which might be due to the Doppler 
modulation effect of the jet and counterjet motion towards 
or away from the observer. 



4. Discussion and conclusions 

As stated in the previous section, if 3C66B really hosts 
an SMBHB in its center, a periodic behavior of the light 
curve in the X-ray and/or 7-ray bands and a certain flux 
ratio / are expected. A measure of P bs and / would allow 
getting both the jet Lorentz factor 7l and the system orbit 
eccentricity e through Fig. [3 and Fig. 0] 

However, since the curves in Fig. 0] do not depend on 
the mass ratio q, also measuring / does not make it pos- 
sible to obtain a complete definition of the orbital param- 
eters of the SMBHB. Thus, a third observation, i.e. the 
gravitational wave emission, is necessary to extract the 
value of q. 

First, note that the emission of gravitational radiation 
makes the SMBH B orb it to sh rink w ithin the merging 
time-scale (Peters GUI Fitchett fl987h 



T 



1 1 1 121e 2 A 

167rc 5 f e ° V 1 + ~W ) 



1181/2299 



Mm(M + m) 



,29/10 



(1 - e 2 ) 3 / 2 



<ie, (13) 



where r is a constant that depends on the initial orbit 
semi- major axis and eccentricity eg. 

In Fig. [5] the coalescing time scale T for different val- 
ues of the mass ratio q (different lines) and two values of 
the semi-major axis a (different colors) is shown as a func- 
tion of the orbit eccentricity e. Models corresponding to 
the higher value of a (red lines) are clearly unlikely since 
the corresponding coalescing time-scale is shorter than 
10 2 yrs. On the contrary, models with the smaller semi- 
major axis (corresponding to that derived by using the 
8.4 GHz map) and q < 10™ 3 always have coalescing time- 
scales longer than ~ 10 2 yrs and may therefore be con- 
sidered acceptable. Note also that, as previously stated, 
since future radio observations at frequencies higher than 
8.4 GHz will allow better constraint of the semi-major 
axis of the SMBHB system the true coalescing time-scale 
of the massive black hole binary systems may be longer 
than that quoted here. 

The SMBHB mass ratio q may be determined if one 
could observe the gravitational wave spectrum. In Fig. 
we present the expected gravitational wave spectrum for 
different values of the mass ratio q and orbital eccentricity 
e. Straight lines correspond to the detection thresholds of 
the next g enera tion of space-based interfer ometer s LISA 
(Rein hard l2000|) and ASTROD (Wu et al. l2000l and Ni 
l2002h . 

Thus, future X-ray and/or 7-ray observations may al- 
low to obtain information about both the bulk motion of 
the jet (tl) and its precession and, consequently, confirm 
the SMBHB model for the radio galaxy 3C66B. These 
measurements may also allow to define the semi-major 
axis a and estimate the orbital eccentricity e of the bi- 
nary system as follows from Figs. |21 and Fig. @] Moreover, 
the mass ratio q may be attained from the observation 
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of the emitted gravitational wave spectrum by using next 
generation space-based interferometers. 

Note that recently a different method for constraining 
the properti es of th e SMBHB in 3C66B has been proposed 
(Jenet et al. l2003lh The main idea is to search for gravi- 
tational wave emission from the SMBHB system by using 
available timing data (7 yrs up t o now ) from the radio 
pulsar PSR B1855+09 (Kaspi et al. ll994|h In fact, gravita- 
tional waves will induce periodic oscillations in the arrival 
times of the individual pulses of the pulsar. Numerical sim- 
ulations performed for different parameters have shown 
that the proposed method can be used to solve the pa- 
rameter degeneracy of the SMBHB in 3C66B only for ex- 
tremely small orbital eccentricity values (e < 0.03). In the 
other cases the method can not gi ve a d efinitive answer 
using the current data (Jenet et al. [2£ 



Fig. 6. The coalescing time scale T as a function of the 
true orbit eccentricity e. Black and red lines correspond to 
the orbital maximum semi-major axis a ~ 5.4 x 10 16 (1 + g) 
cm and a ~ 2.9 x 10 17 (l+g) cm, respectively. Solid, dotted, 
dashed, dot-dashed and two-dot-dashed lines are for q = 
1Q- 5 , 10~ 4 , 10" 3 , 10~ 2 , 10-\ respectively. 




Fig. 7. The expected gravitational wave spectrum for 
the SMBHB in 3C66B. Black, blue and red lines are 
for q = 10~ 5 ,10~ 4 , 10~ 3 , respectively. Solid and dot- 
ted lines indicate eccentricity values of e = 0.7 and 
e = 0.9. Here we have assumed the semi-major axis to 
be a ~ 5.4 x 10 16 (1 + q) cm. Solid, dashed and dotted 
straight lines show the planned sensitivity thresholds of 
LISA, ASTROD1 and ASTROD2 instruments. Here we 
have assumed an integration time of 5 yrs. 
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